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Abstract 

The major objective of this study is to develop an in- 
tegrated orbit, attitude, and structural control system 
architecture for very large Space Solar Power Satellites 
(SSPS) in geosynchronous orbit. This study focuses on 
the 1.2-GW “Abacus” SSPS concept characterized by a 
3.2 x 3.2 km solar-array platform, a 500-m diameter mi- 
crowave beam transmitting antenna, and a 500 x 700 m 
earth-tracking reflector. For this baseline Abacus SSPS 
configuration, we derive and analyze a complete set of 
mathematical models, including external disturbances 
such as solar radiation pressure, microwave radiation, 
gravity-gradient torque, and other orbit perturbation ef- 
fects. The proposed control system architecture utilizes 
a minimum of 500 1-N electric thrusters to counter, si- 
multaneously, the cyclic pitch gravity-gradient torque, 
the secular roll torque caused by an offset of the center- 
of-mass and center-of-pressure, the cyclic roll/yaw mi- 
crowave radiation torque, and the solar radiation pres- 
sure force whose average value is about 60 N. 

1 Introduction 

A renewed interest in space solar power is spurring 
a reexamination of the prospects for generating large 
amounts of electricity from large-scale, space-based so- 
lar power systems. Peter Glaser (Refs. 1-2) first pro- 
posed the Satellite Solar Power Station (SSPS) concept 
in 1968 and received a U.S. patent on a conceptual de- 
sign for such a satellite in 1973. As a result of a series 
of technical and economic feasibility studies by NASA 
and Department of Energy in the 1970s, an SSPS ref- 
erence system was developed in the late 1970s. The 
1979 SSPS reference system, as it is called, featured a 
very large solar array platform (5.3 x 10.7 km) and a 
double-gimballed microwave beam transmitting antenna 
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(1 km diameter). The total mass was estimated to be 
50 x 10 6 kg. A ground or ocean-based rectenna (recti- 
fying antenna) measuring 10 x 13 km would receive the 
microwave beam on the earth and deliver up to 5 GW 
of electricity. 

In 1995, NASA revisited the Space Solar Power (SSP) 
concept to assess whether SSP-related technologies had 
advanced enough to alter significantly the outlook on the 
economic and technical feasibility of space solar power. 
The “Fresh Look” study (Ref. 3), conducted by NASA 
during 1995-1997, found that in fact a great deal had 
changed and that multi-megawatt SSP satellites appear 
viable, with strong space applications. The study also 
found that ambitious research, technology development 
and validation over a period of perhaps 15-20 years 
are required to enable SSP concepts to be considered 
“ready” for commercial development. 

Recent studies by NASA as part of the SSP Ex- 
ploratory Research and Technology (SERT) program 
have produced a variety of new configurations of Space 
Solar Power Satellites (SSPS), including the “Abacus” 
configuration, as described in Refs. 4-6. Some of these 
configurations, such as the “Sun Tower” configuration, 
are based on the passive gravity-gradient stabilization 
concept. However, most other configurations require 
three-axis attitude control to maintain continuous sun 
tracking of the solar arrays in the presence of external 
disturbances including the gravity-gradient torque. A 
cylindrical configuration, which is not affected by the 
troublesome pitch gravity-gradient torque, has also been 
considered by NASA (Ref. 6). 

This study focuses on the 1.2-GW “Abacus” satellite 
configuration shown in Figure 1. This Abacus satellite is 
characterized by its simple configuration consisting of an 
inertially oriented, 3.2 x 3.2 km solar-array platform, a 
500-m diameter microwave beam transmitting antenna 
fixed to the platform, and a 500 x 700 m rotating re- 
flector that tracks the earth. Some unique features of 
the Abacus satellite relative to the 1979 SSPS reference 
system are: 

• The transmitting antenna is not gimballed; instead, 
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Figure 1: Baseline 1.2-GW “Abacus” satellite configu- 
ration (Refs. 4-6). 

an azimuth roll-ring mounted, rotating reflector 
provides earth pointing of the microwave beam. 

• The rotating reflector design thus eliminates mas- 
sive rotary joint and slip rings of the 1979 SSPS 
reference system. 

• Links activated by ball-screw mechanisms tilt the 
reflector to point to ground stations at various lat- 
itudes 

The objectives of this study, in support of the SERT 
program of NASA, are: (i) to develop preliminary con- 
cepts for orbit, attitude, and structural control of very 
large SSPS using a variety of actuators such as control 
moment gyros, momentum wheels, and electric propul- 
sion thrusters; (ii) to develop mathematical models, de- 
fine a top-level control system architecture, and perform 
control system design and analysis for a baseline Aba- 
cus satellite configuration in geosynchronous orbit; and 
(iii) to determine the required number, size, placement, 
mass, and power for the actuators to control the orbit, 
attitude, and structural motions of the baseline Abacus 
satellite. 

2 Mass Properties, Disturbances, 
and Control Requirements 

2.1 Geometric Properties 

The three major parts of the Abacus satellite and their 
dimensions are shown in Figure 1; the mass of each part 


Table 1: Geometric and mass properties of the 1.2-GW 
Abacus satellite 


Solar array mass 

21 xlO 6 kg 

Transmitting antenna mass 

3 x 10 6 kg 

Reflector mass 

0.8 x 10 6 kg 

Total mass 

m = 25 x 10 6 kg 

Platform area 

A = 3200 m x 3200 m 

Area-to-mass ratio 

A/m = 0.4 m 2 /kg 

Roll inertia 

Ji = 2.8 x 10 13 kg-m 2 

Pitch inertia 

J 2 = 1.8 x 10 13 kg-m 2 

Yaw inertia 

J 3 = 4.6 x 10 13 kg-m 2 

cm-cp offset 

200 m (along pitch axis) 

cm-cp offset (uncertainty) 

±20 m (along roll axis) 

Table 2: Solar pressure and microwave radiation distur- 

bances 


Solar pressure force 

(4.5E-6)(1.3)(A) = 60 N 

Solar pressure torque (roll) 

60 N x 200 m 

Solar pressure torque (pitch) 

60 N x 20 m 

Reflector radiation force 

7 N (rotating force) 

Reflector radiation torque 

7 N xl700 m 


is given in Table 1, together with the total mass and area 
of the spacecraft. The mass of the reflector is approx- 
imately 3% of the total mass; therefore, the reflector’s 
mass can be neglected in the analysis of attitude motion, 
simplifying the task in two important respects. First, the 
Abacus satellite can be treated as a single body rather 
than a multibody spacecraft. When the Abacus satellite 
is regarded as rigid, the spacecraft’s moments and prod- 
ucts of inertia for a set of axes fixed in the solar array 
do not vary with time. Second, when the unsymmetrical 
mass distribution of the reflector is left out of account, 
the principal axes of inertia of the spacecraft with re- 
spect to the spacecraft’s mass center are parallel to the 
roll, pitch, and yaw axes illustrated in Figure 1. The 
moments of inertia for these axes, henceforth considered 
to be principal moments of inertia, are given in Table 
1. The center of pressure is located 100 m below the 
geometric center of the square platform, the center of 
mass is located 300 m below the geometric center along 
the pitch axis, and ±20% overall uncertainty in the mass 
properties should be considered in control design. 

2.2 External Disturbances 

External disturbances acting on the Abacus satellite in- 
clude: solar radiation pressure force, microwave radia- 
tion force, gravity-gradient torque, and other orbit per- 
turbation forces. Some of these disturbances with ±20% 
overall uncertainties are summarized in Table 2. Distur- 
bance torques in units of N-m, due to solar pressure, 
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Table 3: Orbit parameters and control requirements 


Table 4: A large single-gimbal CMG 




Cost 

$1M 

Earth’s gravitational parameter 

H = 398,601 km 3 /s 2 

Momentum 

7,000 N-m-s 

Geosynchronous orbit (e, * ~ 0) 

a = 42,164 km 

Max torque 

4,000 N-m 

Orbit period 

23 hr 56 min 4 sec 

Peak power 

500 W 

Orbit rate 

n = 7.292 xlO -5 rad/sec 

Mass 

250 kg 

Stationkeeping accuracy 

±0.1 deg 

Momentum/mass 

28 N-m-s/kg 

Solar array pointing accuracy 

±0.5 deg for roll/pitch 



Microwave beam pointing accuracy 

±5 arcmin 




microwave radiation, cm-cp offset, and cm-cp offset un- 
certainty, can be expressed along the platform-fixed con- 
trol axes as: 


counter the cyclic gravity-gradient torque simply be- 
comes 

3n 2 

u 2 = — ( J3 — J\) sin 2 nt (3) 


Roll: 

di s 

s 12,000 - 11, 900 cos nt 

(la) 

Pitch: 

d 2 s 

s 1200 

(lb) 

Yaw: 

d 3 s 

s — 11, 900 sin nt 

(lc) 


where n is the orbital rate of the Abacus satellite and t 
is time. The constant pitch disturbance torque of 1200 
N-m is due to the assumed cm-cp offset of 20 m along 
the roll axis, and ±20% uncertainty in this disturbance 
model should also be considered in control design. In 
addition to these disturbances, gravity-gradient distur- 
bance torques are also acting on the Abacus satellite. 
It is assumed that the electric currents circulate in the 
solar array structure in such a way that magnetic fields 
cancel out and the Abacus satellite is not affected by the 
magnetic field of the earth. 

2.3 Orbit Parameters and Control Re- 
quirements 

Basic orbital characteristics and control requirements for 
the Abacus satellite in geosynchronous orbit are summa- 
rized in Table 3. 

3 Technical Issues 

3.1 Momentum Storage Requirement 

Assuming that the gravity gradient torque is the only 
external disturbance torque acting along the pitch axis, 
we consider the pitch equation of motion of a rigid space- 
craft in geosynchronous orbit given by 

3n 2 

J2&2 = ~^(J 3 ~ Ji) S1112O2 + U2 (2) 

where Ji, J2, and J3 are, respectively, the roll, pitch, 
and yaw principal moments of inertia; 9 2 is the pitch 
angle measured from the LVLH (local vertical and local 
horizontal) reference frame; n is the orbit rate; and n 2 
is the pitch control torque. 

For continuous sun pointing of the Abacus platform 
with 62 = nt , the pitch control torque required to 


with peak values of ±143,000 N-m. If angular mo- 
mentum exchange devices, such as momentum wheels 
(MWs) or control moment gyros (CMGs), are to be em- 
ployed for pitch control, the peak angular momentum to 
be stored can then be estimated as 

3 n 

ffmax = y (J 3 - Ji) = 2 x 10 9 N-m-s (4) 

This is is about 100,000 times the angular momentum 
storage requirement of the International Space Station 
(ISS). The ISS is controlled by four double-gimballed 
CMGs with a total momentum storage capability of 
about 20,000 N-m-s. The double-gimballed CMGs em- 
ployed by the ISS have a momentum density of 17.5 N- 
m-s/kg, and future advanced flywheels may have a larger 
momentum density of 150 N-m-s/kg. Basic characteris- 
tics of a large single-gimbal CMG are also summarized 
in Table 4. 

Based on the preceding discussion, it can be concluded 
that a traditional momentum management approach us- 
ing conventional CMGs (or even employing future ad- 
vanced flywheels) is not a viable option for controlling 
very large Space Solar Power Satellites. 

To meet the momentum storage requirement of very 
large SSPS, a concept of constructing large-diameter mo- 
mentum wheels in space has been studied in the late 
1970s (Ref. 7). In an attempt to resolve the angular 
momentum storage problem of large sun-pointing space- 
craft, a quasi-inertial sun-pointing, pitch control concept 
was also developed by Elrod in 1972 (Ref. 8), and fur- 
ther investigated by Juang and Wang in 1982 (Ref. 9). 
However, such a “free-drift” concept is not a viable op- 
tion for the Abacus satellite because of the large pitch 
attitude peak error of 18.8 deg and its inherent sensi- 
tivity with respect to initial phasing and other orbital 
perturbations. 

Because the pitch gravity-gradient torque becomes 
naturally zero for cylindrical, spherical or beam-like 
satellites with J\ = J 3 , a cylindrical SSPS configura- 
tion was also studied by NASA to simply avoid such a 
troublesome pitch gravity-gradient torque problem. 
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Figure 2: Solar radiation pressure force acting on an 
ideal flat surface. 


3.2 Solar Radiation Pressure and Large 
Area-to-Mass Ratio 

Despite the importance of the cyclic pitch gravity- 
gradient torque, this study shows that the solar radiation 
pressure force is considerably more detrimental to con- 
trol of the Abacus satellite (and also other large SSPS) 
because of an area-to-mass ratio that is very large com- 
pared to contemporary, higher-density spacecraft. 

The significant orbit perturbation effect of the solar 
radiation pressure on large spacecraft with large area-to- 
mass ratios has been investigated by many researchers 
in the past. A detailed physical description of the solar 
radiation pressure can be found in a recent book on solar 
sailing by Mclnnes (Ref. 10). 

The solar radiation pressure forces are due to photons 
impinging on a surface in space, as illustrated in Figure 
2. Assuming that a fraction, p a , of the impinging pho- 
tons is specularly reflected, a fraction, pd , is diffusely 
reflected, and a fraction, p a , is absorbed by the surface, 
we have 

Ps + Pd + Pa = 1 (5) 

The solar radiation pressure (SRP) force acting on an 
ideal flat surface is then expressed as 

2 p s (n ■ s) + raj 
(6) 

where P = 4.5 x 10 6 N/m 2 is the nominal solar ra- 
diation pressure constant, A is the surface area, n is a 
unit vector normal to the surface, and s is a unit vector 
pointing from the sun to satellite. 

For an ideal case of a perfect mirror with Pd = Pa = 0 
and p s = 1, we have 

F = 2PAcos 2 (f> n 

where Acoscj) is called the projected area of the surface 
under consideration. Also for an ideal case of a black 
body with p s = Pd = 0 and p a = 1, we have 

F = PA cos <j> s 


F = PA(n ■ s) < (p a + p d )s + 


For most practical cases of satellites with small angles 
of </>, the SRP perturbation force per unit mass is simply 
modeled as 

f = P(l + p)(A/m)s (7) 

where p is the overall surface reflectance (0 for a black 
body and 1 for a mirror) and A/m is the area-to-mass 
ratio. 

Let the SRP perturbation acceleration be expressed 
as 

/= f r e r + fees + f z e z (8) 

where {e r , eg, e z } is a set of unit vectors of the so-called 
perifocal reference frame. Ignoring the effects of sea- 
sonal variations of the sun vector, we simply obtain 
f r « fsind and fg ~ f cos6 where / = P(1 + p)(A/m ) 
and 9 is the true anomaly. 

From the orbit perturbation analysis (Refs. 11-12), 
we have 


da 

dt 

de 

dt 


n\J 1 — e 
Vl - i 


z{f r esm6 + fg(l + ecos0)} (9a) 


na 


~{f r sin 9 + f e (cos9 + cosE)} (9b) 


where 9 and E are the true and eccentric anomalies, 
respectively. For geosynchronous satellites with e « 0, 
we obtain 


and 


da 2 2/ . 

-rr = ~fe = — sin ( 
at n n 

=> A a = 0 per day 


dc 

— = — (f r sin 9 + 2/ e cos 9) 
at na 

= — (/ sin 2 9 + 2/ cos 2 9) 
na 

= — ( x + x cos 29 


Ae ; 


( 10 ) 


3t r/ 

per day 

ma 


( 11 ) 


The solar radiation pressure effect on the longitude 
change can also be found as 

3n 2 


^ dn 3 n da 

dt 2 a dt 

= - 3/ cos , 


,„=- 3 f e 

Za n a 


(12) 


For the Abacus satellite, we have 


Area-to-mass ratio A/m ss 0.4 m 2 /kg 
SRP force « 60 N 

SRP perturbation acceleration « 2.4 x 10^ 6 
3 tt( 4.5 x 10~ 6 )(1.3)A/m 


Ae 


I I i i 


/ ^ 
m/s 

~ 4 per day 


Longitude drift AA = 2Ae « 0.0115 deg/day 

maximum Ae ss 0.018 at the mid year 

maximum AA = 2Ae « 2 deg; maximum Aa « 1.8 km 
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Figure 4: A large Space Solar Power Satellite (SSPS) in 
geosynchronous orbit. 


Figure 3: Orbit simulation results of the Abacus satellite 
with the effects of the earth’s oblateness and triaxiality, 
luni-solar perturbations, and 60-N solar radiation pres- 
sure force (time in units of orbits) . 


Consequently, orbit eccentricity control using high-/ sp 
ion engines becomes necessary. The yearly propellant re- 
quirement to counter the approximately 60-N solar pres- 
sure force can be estimated as 


A (60) (24 x 3600 x 365) 

m ~ 5000 x 9.8 


40, 000 kg/year 


correspond to an initial position and velocity of 

f= XI + YJ+ ZK = -11698.237/ - 40508.869 J (km) 
v = 2.954/ — 0.853 J (km/s) 

where {/, J, K} is a set of unit vectors of the Earth- 
Centered-Inertial (ECI) reference frame. 

In the next section, we develop an attitude dynam- 
ics model of sun-pointing spacecraft in geosynchronous 
orbit for attitude control system architecture design. 


where an ion engine with I sp of 5000 sec is assumed. De- 
tailed discussions of an electric propulsion system pro- 
posed for the Abacus satellite will be presented later in 
this paper. 

Typical north-south and east- west stationkeeping ma- 
neuvers for the Abacus satellite will also require 

Am = to ^1 — exp ~ 30,000 kg/year 

where to = 25 x 10 6 kg, AV = 50 m/s per year, g = 9.8 
m/s 2 , and I sp = 5000 sec. 

The results of 30-day simulations of orbital motion 
of the Abacus satellite, with the effects of the earth’s 
oblateness and triaxiality, luni-solar perturbations, and 
60-N solar pressure force included, are shown in Figure 
3. It is worth noting the extent to which eccentricity 
and inclination are perturbed. 

The initial values used in the simulations corre- 
spond to a circular, equatorial orbit of radius 42164.169 
km; therefore, the initial orbital elements are: a = 

42164.169 km and e = i = H= u>~0. The epoch used 
to calculate the solar and lunar positions, as well as the 
earth’s orientation in inertial space, is March 21, 2000. 
In order to place the spacecraft at an initial terrestrial 
longitude of 75.07 deg (one of the stable longitudes), a 
true anomaly 9 of 253.89 deg is used. These elements 


4 Attitude Equations of Motion 
of Sun-Pointing Spacecraft 

Consider a spacecraft in circular orbit, as illustrated in 
Figure 4. A local vertical and local horizontal (LVLH) 
reference frame A with its origin at the center of mass 
of an orbiting spacecraft has a set of unit vectors 
{ 01 , 02 , 03 } with oi along the orbit direction, 02 per- 
pendicular to the orbit plane, and 03 toward the earth, 
as illustrated in Figure 4. The angular velocity of A with 
respect to an inertial or Newtonian reference frame N is 

lu a / n = —na2 (13) 

where n is the constant orbital rate. The angular ve- 
locity of the body-fixed reference frame B with basis 
vectors {&i, 62 , ^ 3 } is then given by 

jB/N = jB/A + jA/N = Q B/A _ ^ (14) 

where Q B I A is the angular velocity of B relative to A. 

To describe the orientation of the body-fixed reference 
frame B with respect to the LVLH reference frame A in 
terms of three Euler angles Qi ( i = 1,2,3), consider the 
sequence of Ci(0i) <— 03 (^ 3 ) <— 02 (^ 2 ) from the LVLH 
reference frame A to a body- fixed reference frame B. 
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For this rotational sequence, we have 


bi 


' C n 

C12 

c 13 ' 


Si 


i 2 

= 

C21 

C22 

C23 


a 2 

( 15 ) 

b 3 


C31 

C32 

C33 


. . 



J2W2 — ( J 3 — Ji)lo 3 loi — —3n 2 (J 3 — J\)C 33 Ci 3 
J&3 — ( J\ — J2)u>ii02 = —3 n 2 (Ji — J 2 )Ci 3 C 2 3 

( 20 ) 

where 


where C\\ = cos 02 cos $3, C12 = sin 0 3 , C13 = 
— sin 9 2 cos 03, etc. 

Assuming that the gravity gradient torque is the only 
external disturbance torque acting along the pitch axis, 
we obtain the rotational equation of motion of a rigid 
body with an angular momentum H in a circular orbit 
as 

= “ <“> 

where {d/dt}]y indicates differentiation with respect to 
time in reference frame N and {d/dt}s indicates differ- 
entiation with respect to time in reference frame B. The 
gravity-gradient torque M is expressed in vector-dyadic 
form as: 

M = 3n 2 a 3 x J ■ CL3 ( 17 ) 

where n = \J is the orbital rate, 03 = —R C /R C1 and 
J is the inertia dyadic of the spacecraft with respect to 
its mass center (Refs. 12 13 ). 

Since H = J -Co B / N , the attitude dynamical equations 
of motion can be rewritten as 

J ■ Co + Co x J ■ Co = 3n 2 a 3 x J ■ a 3 ( 18 ) 

where Co = Co B//N and Co = { dCo/dt = {dCo/dt}s- Ex- 
pressing Co, a 3 , and J in terms of basis vectors of the 
body-fixed reference frame B as 


Co = ooibi + W262 + oo 3 b 3 

( 19 a) 

U3 = C1361 + C23&2 + C3363 

( 19 b) 

3 3 

j = yy JijbSj 

( 19 c) 


*= l j=i 


we obtain the attitude dynamical equations of motion in 
matrix form as: 






Jll 

J12 

d\ 3 

C01 





J21 

J22 

J23 

C02 





J31 

J32 

J33 

. 


0 

-<03 

L0 2 

Jll 

J12 

d[ 3 

COi 

+ 

lo 3 

0 

-LOl 

J21 

J22 

J23 

002 


—L02 

LOi 

0 

J31 

J32 

J33 

. W 3 



0 

~c 33 

C23 

Jll 

J12 

Jl3 

Cl 3 

3 n 2 

C33 

0 

-c 13 

J21 

J22 

J23 

C23 


1 

CO 

C13 

0 

J31 

J32 

J33 

C33 


The dynamical equations of motion about the body- 
fixed principal axes become 

J\C0\ — (J2 — J 3 )i02(O 3 = — 3n 2 (J2 — J 3 )C 23 C 33 


Ci 3 = — sin 0 2 cos 0 3 

C23 = cos 6\ sin 0 2 sin 0 3 + sin 9\ cos 9 2 

C33 = — sin 9\ sin 9 2 sin 9 3 + cos 9 1 cos 02 

for the sequence of C\(0\) <— C 3 (0 3 ) <— C 2 ($2) under 
consideration. For this rotational sequence we have the 
following kinematic differential equations: 


C 0 3 

— C 0! S 0 3 

S 0! S 0 3 

COi 


' 0 ' 

0 

C 0i 

— S0i 

L0 2 

+ 

n 

0 

S 01 C 03 

C 0i C 0 3 

. W 3 


0 


(21) 


where c 9i = cos 9i and s 9i= sin 9i . 

One may linearize Eqs. ( 20 )-( 21 ) “about” an LVLH 
orientation while admitting a large pitch angle as follows. 
Assume 9\ and 9 3 remain small, allow 0 2 to be large, 
assume uo\ and co 3 are small, and L02 is equal to the sum 
of a small quantity and —n. The attitude equations of 
motion that are linear in the small quantities can then 
be obtained as (Ref. 12 ): 

J\9\ + (1 + 3cos 2 0 2 )u 2 (J 2 — J^)9\ — n( Ji — J 2 + J%)9 3 
+ 3( J 2 — J 3 )n 2 (sin 9 2 cos 9 2 )9 3 = u\ + di 
J202 + 3n 2 ( Ji - J 3 ) sin 0 2 cos 9 2 = u 2 + d 2 
J 3 9 3 + (1 + 3sin“02)u 2 (J2 — Ji)9 3 + n(J\ — J 2 + J 3 )9± 
+ 3 ( J 2 - Ji)n 2 (sin 9 2 cos 9 2 )9 1 = u 3 + d 3 

( 22 ) 

where u r and d, are control and disturbance torques, 
respectively. 

However, for a quasi-inertially stabilized, sun-pointing 
SSPS in geosynchronous orbit with small body rates, u>i 
(i = 1 , 2 , 3 ), and small roll/yaw angles, 9\ and 9 3 , the 
kinematic differential equations, (21), can be linearized 


in the small quantities, as 

follows: 


0i ^ 

! <Ol 

( 23 a) 

02~ 

'■ lo 2 + n 

( 23 b) 

03 ~ 

: LO 3 

( 23 c) 


The attitude equations of motion of a quasi-inertially 
stabilized, sun-pointing spacecraft with small roll and 
yaw angles, 0\ and 9 3 , can then be found as 


J\9\ — 

-3 n 2 (J 2 

- J 3 )(cos 2 0 2 )0i 


- 3 ( J 2 - 

J 3 )n 2 ( sin 0 2 cos 9 2 )9 3 + u\ + d\ 

•^202 = 

- 3 n 2 (Ji 

- J 3 ) sin 0 2 cos 0 2 + u 2 + d 2 

J363 = 

-3 n 2 (J 2 

- Ji)(sin 2 0 2 )0 3 


-3 (J 2 - 

Ji)n 2 (sin 0 2 cos 0 2 )0i + ^3 + ^3 


( 24 ) 
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Node Number Locations for Normal Modes Results 

(Nine Nodes Shown in Red) 

4 3200m (80 arrays) 



Transmitter t X Reflector (500m x 750m) 

(500m diameter) 

Front View (Abaqus support truss in back) 






Figure 6: Bode magnitude plots of reduced-order trans- 
fer functions from an input force at node #1 to various 
output locations. 


Figure 5: Selected FEM node locations for control analy- 
sis and design (Courtesy of Tim Collins at NASA LaR.C). 

The pitch attitude angle relative to the LVLH frame, 02, 
is not restricted to be small, but it may be regarded as 
a sum, 02 = rot + <502, where <5 0 2 is a small pitch attitude 
error. Kinematical and dynamical differential equations 
can then be made linear in the small quantities uq, ui 2 , 
co 3 , 0i, <50 2 , and 0 3 . For such a case, Eqs. (21) become 

0i « u q; <50 2 « w 2 ; 03 ~ w 3 (25) 

and Eqs. (24) become 

Ji0i + 3ro 2 (J2 — J 3 )[(cos 2 rof)0i + (l/2)(sin 2rot)0 3 ] 

= ui + di 

J 2 <502 + 3ro 2 (Ji — J 3 )[(cos 2 rot — sin 2 rot)<50 2 
+(1/2) sin2rot] = u 2 + d 2 
J 3 0 3 + 3ro 2 (J2 — +i)[(sin 2 rot)0 3 + (l/2)(sin 2rot)0i] 

= u 3 + d 3 

(26) 

where (0 i,< 502,0 3 ) are the small roll, pitch, and yaw at- 
titude errors of a sun-pointing spacecraft, respectively. 

Equations (24) or (26) are the attitude equations of 
motion of the Abacus satellite for control design in the 
presence of the external disturbances, di, described as in 
Eqs. (1). 

5 Structural Dynamic Models 

Dynamics and control problems of large flexible plat- 
forms in space, such as the square Abacus satellite, have 
been investigated by many researchers in the past. The 
flexible structural dynamics and control problem is a 


topic of continuing practical as well as theoretical inter- 
est. However, a significant control-structure interaction 
problem, possible for such a very large Abacus platform 
(3.2 x 3.2 km) with the lowest structural mode frequency 
of about 0.002 Hz, is avoided simply by designing an at- 
titude control system with very low bandwidth (< orbit 
frequency of 1 x 10 -5 Hz). The proposed low-bandwidth 
attitude control system (to be presented in the next sec- 
tion), however, utilizes a concept of cyclic-disturbance 
accommodation control to provide the required ±5 ar- 
cmin pointing accuracy of the Abacus platform in the 
presence of large, but slowly varying, external distur- 
bances and dynamic modeling uncertainties. Conse- 
quently, the flexible structural control problem is not 
further elaborated in this study, while a structural dy- 
namic interaction problem with thermal distortion needs 
to be further investigated in a future study. 

Various structural concepts for providing the required 
stiffness and rigidity of the Abacus platform have been 
presented in Refs 4-5. Selected node locations for con- 
trol analysis and design are shown in Figure 5. Typ- 
ical pole-zero patterns of reduced-order transfer func- 
tions can be seen in Figure 6. Computer simulation re- 
sults of a reduced-order structural model with the lowest 
16 modes, confirm that the control-structure interaction 
problem can be simply avoided by a low-bandwidth at- 
titude control system. 

6 Control System Architecture 

The area-to-mass ratio of the Abacus satellite, A/m = 
0.4 m 2 /kg, relatively large when compared to 0.02 m 2 /kg 
of typical geosynchronous communications satellites, is 
a key parameter characterizing the very large size of the 
Abacus satellite. If left uncontrolled, this can cause a 
cyclic drift in the longitude of the Abacus satellite of 2 
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deg, east and west. Thus, in addition to standard north- 
south and east- west stationkeeping maneuvers for ±0.1 
deg orbit position control, active control of the orbit ec- 
centricity using ion thrusters with high specific impulse, 
I sp , becomes mandatory. Furthermore, continuous sun 
tracking of the Abacus satellite requires large control 
torques to counter various disturbance torques. A con- 
trol system architecture developed in this study utilizes 
properly distributed ion thrusters to counter, simultane- 
ously, the cyclic pitch gravity-gradient torque, and solar 
radiation pressure. 

6.1 Electric Propulsion System 

Basic characteristics of an electric propulsion system for 
the Abacus satellite are summarized in Table 5. Approx- 
imately 85,000 kg of propellant per year is required for si- 
multaneous orbit, attitude, and structural control using 
a minimum of 500 1-N electric propulsion thrusters with 
I sp = 5,000 sec. The yearly propellant requirement is re- 
duced to 21,000 kg if an I sp of 20,000 sec can be achieved 
(as was assumed for the 1979 SSPS reference system). 
As I sp is increased, the propellant mass decreases but 
the electric power requirement increases; consequently, 
the mass of solar arrays and power processing units in- 
creases. Based on a minimum of 500 1-N thrusters, a 
mass/power ratio of 5 kg/kW, and a power/thrust ra- 
tio of 30 kW/N, the total dry mass (power processing 
units, thrusters, tanks, feed systems, etc.) of an electric 
propulsion system proposed for the Abacus satellite is 
estimated as 75,000 kg. 

The capability of present electric thrusters is orders 
of magnitude below that required for the Abacus satel- 
lite. If the xenon fueled, 1-kW level, off-the-shelf ion en- 
gines available today are to be employed, the number of 
thrusters would be increased to 15,000. The actual total 
number of ion engines will further increase significantly 
when we consider the ion engine’s lifetime, reliability, 
duty cycle, and redundancy. 

For example, the 2.3-kW, 30-cm diameter ion engine 
of the Deep Space 1 spacecraft has a maximum thrust 
level of 92 mN. Throttling down is achieved by feeding 
less electricity and xenon propellant into the propulsion 
system. Specific impulse ranges from 1,900 sec at the 
minimum throttle level to 3,200 sec. 

In principle, an electric propulsion system employs 
electrical energy to accelerate ionized particles to ex- 
tremely high velocities, giving a large total impulse for 
a small consumption of propellant. In contrast to stan- 
dard propulsion, in which the products of chemical com- 
bustion are expelled from a rocket engine, ion propulsion 
is accomplished by giving a gas, such as xenon (which 
is like neon or helium, but heavier), an electrical charge 
and electrically accelerating the ionized gas to a speed 
of about 30 km/s. When xenon ions are emitted at such 
high speed as exhaust from a spacecraft, they push the 
spacecraft in the opposite direction. However, the ex- 


Table 5: Electric propulsion system for the 1.2-GW Aba- 
cus satellite 


Thrust, T 

> 1 N 

Specific impulse, I sp = T/(fng) 

> 5,000 sec 

Exhaust velocity, V e = I sp g 

> 49 km/s 

Total efficiency, g = P 0 /Pi 

> 80% 

Power/thrust ratio, Pi/T 

< 30 kW/N 

Mass/power ratio 

< 5 kg/kW 

Total peak thrust 

200 N 

Total peak power 

6 MW 

Total average thrust 

80 N 

Total average power 

2.4 MW 

Number of 1-N thrusters 

> 500 

Total dry mass 

> 75,000 kg 

Propellant consumption 

85,000 kg/yeai 


Note: T = mV e , P a = \mV? = \TV e , P 0 /T = ^V e = 
ideal power/thrust ratio, Pi/T = ^ l -V e , I sp = T/(mg) = 
V e /g , V e = I sp g where g = 9.8 m/s 2 , to is the exhaust 
mass flow rate, Pi is the input power, and P 0 is the 
output power. 

haust gas from an ion thruster consists of large num- 
bers of positive and negative ions that form an essen- 
tially neutral plasma beam extending for large distances 
in space. It seems that little is known yet about the 
long-term effect of such an extensive plasma on geosyn- 
chronous satellites with regard to communications, solar 
cell degradation, environmental contamination, etc. 

6.2 Control System Architecture 

A control system architecture developed in this study is 
shown in Figure 7. The proposed control system utilizes 
properly distributed ion thrusters to counter, simultane- 
ously, the cyclic pitch gravity-gradient torque, the secu- 
lar roll torque caused by an offset of the center-of-mass 
and center-of-pressure, the cyclic roll/yaw microwave ra- 
diation torque, and the solar radiation pressure force 
whose average value is about 60 N. 

A significant control-structure interaction problem, 
possible for such very large Abacus platform with the 
lowest structural mode frequency of 0.002 Hz, is simply 
avoided by designing an attitude control system with 
very low bandwidth (< orbit frequency). However, the 
proposed low-bandwidth attitude control system utilizes 
a concept of cyclic-disturbance accommodating control 
to provide ±5 arcmin pointing of the Abacus platform 
in the presence of large, but slowly varying, external dis- 
turbances and dynamic modeling uncertainties. 

The concept of cyclic-disturbance accommodation 
control has been successfully applied to a variety of 
space vehicle control problems (Ref. 12), including the 
Hubble Space Telescope, the International Space Sta- 
tion, flexible space structures, and halo orbit control. 
High-bandwidth, colocated direct velocity feedback, ac- 



Figure 7: An integrated orbit, attitude, and structural 
control system architecture employing electric propul- 
sion thrusters. 

tive dampers may need to be properly distributed over 
the platform. Detailed technical discussions of practical 
control systems design for space vehicles in the presence 
structural flexibility as well as persistent external dis- 
turbances can be found in Ref. 12. Thus, theoretical 
aspects of the control law design problem of the Abacus 
satellite are not elaborated upon in this paper. 

Placement of approximately 500 1-N electric propul- 
sion thrusters at 12 different locations is illustrated in 
Figure 8. 

In contrast to a typical placement of thrusters at the 
four corners, e.g., employed for the 1979 SSPS refer- 
ence system, the proposed placement shown in Figure 
8 minimizes roll/pitch thruster couplings as well as the 
excitation of platform out-of-plane structural modes. A 
minimum of 500 ion engines of 1-N thrust level are re- 
quired for simultaneous attitude and stationkeeping con- 
trol. When reliability, lifetime, duty cycle, lower thrust 
level, and redundancy of ion engines are considered, this 
number will increase significantly. 

6.3 Control Simulation Results 

Control system simulation results of an illustrative case 
with 10-deg initial attitude errors are shown in Figure 
9. In this simulation, various dynamic modeling un- 
certainties (i.e., ±20 % uncertainties for the moments 
and products of inertia, cm-cp offset, external distur- 
bances, etc.) are included. The proposed low-bandwidth 


Roll: 1/3 Pitch: 2/4 Yaw: 5/6/7/8 
Orbit Eccentricity, Roll/Pitch Control: 1/3, 2/4 
E/W and Yaw Control: 9/10/11/12 
N/S and Yaw Control: 5/6Z7/8 

Figure 8: Placement of a minimum of 500 1-N electric 
propulsion thrusters at 12 different locations, with 100 
thrusters each at locations #2 and #4. 



Figure 9: Control simulation results for the proposed 
attitude control system architecture. 
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a (km) 



Figure 10: Orbit control simulation results with continu- 
ous (non-impulsive) eccentricity and inclination control 
(time in units of orbits). 

attitude control system, which effectively utilizes the 
concept of cyclic-disturbance accommodation control, 
maintains the required ±5 arcmin, steady-state point- 
ing of the Abacus platform in the presence of large, but 
slowly varying, external disturbances and dynamic mod- 
eling uncertainties. The total thrusting force from the 
roll/pitch thrusters #1 through #4 nearly counters the 
60-N solar radiation pressure force; however, any resid- 
ual AV caused by various uncertainties should be cor- 
rected during standard east-west stationkeeping maneu- 
vers. 

Orbit control simulation results with the effects of 
earth’s oblateness and triaxiality, luni-solar perturba- 
tions, 60-N solar radiation pressure force, and simul- 
taneous orbit and attitude control thruster firings are 
shown in Figure 10. It can be seen that the eccentricity 
and inclination are all properly maintained. The feasi- 
bility of using continuous (non-impulsive) firings of ion 
thrusters for simultaneous attitude and orbit control is 
demonstrated in this study; however, a further detailed 
orbit control study is needed. The orbit control prob- 
lem of geosynchronous satellites is a topic of continuing 
practical interest (Refs. 14-16). 

7 Summary and Recommenda- 
tions 

7.1 Summary of Study Results 

Despite the importance of the cyclic pitch gravity- 
gradient torque, this study has shown that the solar ra- 
diation pressure force is considerably more detrimental 
to control of the Abacus satellite (and also other large 
SSPS) because of an area-to-mass ratio that is very large 


compared to contemporary, higher-density spacecraft. 

A key parameter that characterizes the Abacus satel- 
lite is its area-to-mass ratio, A/m , of 0.4 m 2 /kg, which is 
relatively large when compared to 0.02 m 2 /kg for typical 
geosynchronous communications satellites. 

A significant control-structure interaction problem, 
possible for such very large, flexible Abacus platform 
with the lowest structural mode frequency of 0.002 Hz, 
is simply avoided by designing an attitude control sys- 
tem with very low bandwidth (< orbit frequency). How- 
ever, the proposed low-bandwidth attitude control sys- 
tem effectively utilizes a concept of cyclic-disturbance 
accommodating control to provide ±5 arcmin pointing of 
the Abacus platform in the presence of large, but slowly 
varying, external disturbances and dynamic modeling 
uncertainties. 

Approximately 85,000 kg of propellant per year is re- 
quired for simultaneous orbit, attitude, and structural 
control using a minimum of 500 1-N electric propul- 
sion thrusters with I sp = 5,000 sec. The total dry mass 
(power processing units, thrusters, tanks, feed systems, 
etc.) of an electric propulsion system proposed for the 
Abacus satellite is estimated as 75,000 kg. 

7.2 Recommendations for Future Re- 
search 

The baseline control system architecture developed for 
the Abacus satellite requires a minimum of 500 ion en- 
gines of 1-N thrust level. The capability of present elec- 
tric thrusters is orders of magnitude below that required 
for the Abacus satellite. If the xenon fueled, 1-kW level, 
off-the-shelf ion engines available today, are to be em- 
ployed, the number of thrusters would be increased to 
15,000. The actual total number of ion engines will fur- 
ther increase significantly when we consider the ion en- 
gine’s lifetime, reliability, duty cycle, redundancy, etc. 
Consequently, a 30-kW, 1-N level electric propulsion 
thruster with a specific impulse greater than 5,000 sec 
needs to be developed for the Abacus satellite if an ex- 
cessively large number of thrusters is to be avoided. 

Several high-power electric propulsion systems are 
currently under development. For example, the NASA 
T-220 10-kW Hall thruster recently completed a 1,000- 
hr life test. This high-power (over 5 kW) Hall thruster 
provides 500 mN of thrust at a specific impulse of 2,450 
sec and 59% total efficiency. Dual-mode Hall thrusters, 
which can operate in either high-thrust mode or high- 
I sp mode for efficient propellant usage, are also being 
developed. 

The exhaust gas from an electric propulsion system 
forms an essentially neutral plasma beam extending for 
large distances in space. Because little is known yet 
about the long-term effect of an extensive plasma on 
geosynchronous satellites with regard to communica- 
tions, solar cell degradation, etc, the use of lightweight, 
space-assembled large-diameter momentum wheels may 
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Table 6: Technology advancement needs for the Abacus 
SSPS 


Electric Thrusters 

30 kW, 1 N 
I sp > 5000 sec 
(500-1,000 thrusters) 

CMGs 

2,000 N-m-s/kg 
500,000 N-m-s/unit 

Space- Assembled 

66,000 N-m-s/kg 

Momentum Wheels 

4 xlO 8 N-m-s/unit 
(5-10 MWs) 


also be considered as an option for the Abacus satellite; 
therefore, these devices warrant further study. In Table 
6, the technology advances required for the Abacus satel- 
lite are summarized. It is emphasized that both electric 
propulsion and momentum wheel technologies require 
significant advancement to support the development of 
large SSPS. 

Despite the huge size and low structural frequencies 
of the Abacus satellite, the control-structure interaction 
problem appears to be a tractable one because the tight 
pointing control requirement can be met even with a 
control bandwidth that is much lower than the lowest 
structural frequency. However, further detailed study 
needs to be performed for achieving the required 5- 
arcmin microwave beam pointing accuracy in the pres- 
ence of transmitter/rcflector-coupled structural dynam- 
ics, Abacus platform thermal distortion and vibrations, 
hardware constraints, and other short-term impulsive 
disturbances. 

Although the rotating reflector concept of the Abacus 
satellite eliminates massive rotary joint and slip rings of 
the 1979 SSPS reference concept, the transmitter fixed 
to the Abacus platform results in unnecessarily tight 
pointing requirements imposed on the platform. Further 
system-level tradeoffs will be required for the microwave- 
transmitting antenna design, such as whether or not to 
gimbal it with respect to the platform, use mechanical 
or electronic beam steering, and so forth. 
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